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In an earlier study, an in vitro evolution procedure was applied to a large 
population of variants of the Tetrahymena group I ribozyme to obtain 
individuals with a lCf-fold improved ability to cleave a target 
single-stranded DNA substrate under simulated physiological conditions. 
The evolved ribozymes also showed a twofold improvement, compared to 
the wild-type, in their ability to cleave a single-stranded RNA substrate. 
Here, we report continuation of the in vitro evolution process using a new 
selection strategy to achieve both enhanced DNA and diminished 
RNA-cleavage activity. Our strategy combines a positive selection for 
DNA cleavage with a negative selection against RNA binding. After 36 
"generations" of in vitro evolution, the evolved population showed an 
~10G- fold increase in the ratio of DNA to RNA-cleavage activity. 
Site-directed mutagenesis experiments confirmed the selective advantage 
of two covarying mutations within the catalytic core of the ribozyme that 
are largely responsible for this modified behavior. The population of 
ribozymes has now undergone a total of 63 successive generations of 
evolution, resulting in an average of 28 mutations relative to the wild-type 
that are responsible for the altered phenotype. 

© 1996 Academic Press Limited 

Keywords: catalytic RNA; DNA cleavage; in vitro evolution; ribozyme 


Introduction 

In vitro evolution has proven to be a powerful 
technique for obtaining novel catalytic RNAs (for 
reviews, see Breaker & Joyce, 1994; Lorsch & 
Szostak, 1996). One approach begins with a pool of 
variants of an existing ribozyme that catalyzes a 
reaction similar to the one desired. Repeated 
rounds of selection, amplification, and mutation are 
carried out until the desired catalyst is obtained. 
Previous examples of this approach have demon- 
strated that the original catalytic activity usually is 
retained and may even improve, suggesting that 
in vitro evolution has served to broaden rather than 
transform the activity of the starting enzyme. For 
example, a group I ribozyme that is active in the 
presence of Mg^ 4 but not Ca 2+ was evolved to 
function in the presence of Ca 2 \ but did not lose its 
ability to function in the presence of Mg 2 * (Lehman 
& Joyce, 1993). More pertinent to the present 


Abbreviations used: IGS, internal guide sequence; 
PCR, polymerase chain reaction; Gn, generation after 
n rounds of in vitro evolution; E, ribozyme; rP, RNA 
product; dS, DNA substrate; dP, DNA product; wt, 
wild-type; PI helix, ribozyme/substrate duplex. 


study, a group I ribozyme that efficiently cleaves 
RNA but has a barely detectable ability to cleave 
DNA was evolved to cleave DNA efficiently, yet 
showed a slight improvement in its ability to cleave 
RNA (Beaudry & Joyce, 1992; Tsang & Joyce, 1994). 
Here we describe a novel in vitro evolution scheme 
that simultaneously selects for the desired activity 
and against the existing one, leading to specializ- 
ation rather than expansion of the enzyme's 
capabilities. 

We began this study with variants of the 
self-splicing group I intron of Tetrahymena ther- 
mophila, a ribozyme that catalyzes the sequence- 
specific cleavage of a single-stranded RNA 
substrate via a phosphoester transfer mechanism 
(Figure 1; Zaug & Cech, 1986; Zaug el al. f 1986). The 
ribozyme binds its substrate in two steps: first, a 
region at the 5’ end of the ribozyme, termed the 
"internal guide sequence" (IGS), forms Watson - 
Crick base-pairs with the substrate; then the 
ribozyme/substrate duplex (PI helix) forms ter- 
tiary interactions with residues in the catalytic core 
of the ribozyme. This positions a specific phospho- 
ester bond of the substrate in close proximity to the 
3*-hydroxyl of a bound guanosine nucleoside, 
including the guanosine that lies at the 3’ end of the 
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ribozyme (Herschlag, 1992; Bevilacqua et nL, 1992). 
A phosphoester transfer reaction ensues, involv- 
ing nucleophilic attack by the guanosine 3'-OH, 


and resulting in cleavage of the substrate with 
concomitant ligation of the 3' portion of the 
substrate to the 3 1 oxygen of guanosine (Michel 
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Figure 1. Secondary structure of the Tetrnhymem ribozyme (Cech et nL, 1994). Paired structural elements are 
designated by Pi. Joining regions between paired elements i and /, referred to as J i/j, are not labeled. Nucleotide 
positions that were partially randomized in the initial population are indicated by shaded regions. The internal guide 
sequence (IGS) is shown in bold face, and the DNA substrate is shown in lowercase letters. Nucleotide positions 
discussed in the text are labeled. Positions that were subjected to site-directed mutagenesis are boxed. 
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Figure 2. Selection strategy. The 
ribozyme population (E) was chal- 
lenged to cleave the DNA substrate 
(dS) in the presence of an inhibitor 
corresponding to the 5* product 
portion of the RNA-cleavage reac- 
tion (rP). Ribozymes that cleaved 
the DNA substrate became attached 
to the 3* portion of the substrate, 
which served as a tag for sub- 
sequent selective amplification. Ri- 
bozymes that bound the RNA 
product were prevented from cleav- 
ing the DNA substrate and were not 
amplified. 



et at., 1989; Been & Perotta, 1991; Mei & Herschlag, 
1996). 

Certain 2-hydroxyl groups within the RNA 
substrate form tertiary contacts with residues in 
the catalytic core of the ribozyme during the 
second step of substrate binding (Pyle & Cech, 
1991; Pyle et al., 1992, Bevilacqua & Turner; 1991). 
In addition, the 2'-hydroxyl of the ribonucleo- 
tide at the cleavage site contributes — 10 3 -fold 
to transition-state stabilization during the trans- 
esterification reaction (Herschlag et oL, 1993). It 
is not surprising, therefore, that the wild-type 
ribozyme cleaves a DNA substrate very ineffi- 
ciently, with a ICP-fold reduction in substrate 
binding affinity and a 1 CP-fold reduction in catalytic 
rate compared to the reaction With an analogous 
RNA substrate. 

Beaudry & Joyce (1992) used an in vitro evolution 
procedure to improve the DNA-cleaving activity of 
the Tetrahymena ribozyme by about 100- fold under 
simulated physiological conditions. A pool of 10 13 
variants of the Tetrahymena ribozyme was gener- 
ated, and individuals with superior DNA-cleaving 
activity were preferentially amplified on the basis 
of attachment of the 3' portion of the substrate to 
the 3’ end of the ribozyme, which occurred as a 
consequence of the cleavage reaction. Amplified 
molecules were utilized as templates for cDNA 
synthesis, the resulting cDNA was amplified by a 
polymerase chain reaction (PCR), and the PCR 
products were transcribed to generate a new pool 
of RNAs. This procedure was carried out for nine 
successive generations, resulting in RNA molecules 
that cleaved a single-stranded DNA substrate with 
a k cM of 0.06 min"' and K m of 2 pM. This represents 
an improvement of 30-fold and threefold, respect- 
ively, compared to the wild-type Tetrahymena 
ribozyme. 


In a subsequent in vitro evolution study, Tsang & 
Joyce (1994) set out to optimize both k CM and K m by 
appropriate manipulation of the selection con- 
straints. First, the substrate concentration was 
reduced 50-fold during the target reaction to favor 
ribozymes with high substrate binding affinity. 
This resulted in a 25-fold improvement in both K d 
and K m over nine generations of in vitro evolution. 
Next, the time allowed for the cleavage reaction 
was decreased by 12-fold to favor ribozymes with 
a fast rate of reaction. This resulted in a 50-fold 
improvement in k c * » over nine generations, with no 
significant change in K m . Thus, after a total of 27 
generations of in vitro evolution, ribozymes were 
obtained that cleaved DNA with a catalytic 
efficiency, k CiX /K mt of 3.2 x 10 6 M" 1 min" 1 , represent- 
ing a —10 5 -fold improvement compared to the 
wild-type. However, as noted above, RNA-cleav- 
age activity also improved, so that ribozymes 
obtained after the 27th generation still cleaved RNA 
about tenfold more efficiently than they cleaved 
DNA. Thus, the in vitro evolution process had led 
to a broadening of the catalytic abilities of the 
ribozyme. Nature had evolved an RNA catalyst 
with the ability to cleave RNA; then, in the 
laboratory, this activity was generalized to include 
both RNA and DNA substrates. 

In this study, we devised a novel in vitro 
evolution strategy in order to improve the 
specificity of the ribozyme for DNA compared to 
RNA substrates. Our scheme combined a positive 
selection for DNA cleavage with a negative 
selection against RNA binding. As shown in Fig- 
ure 2, ribozymes were challenged to cleave a DNA 
substrate in the presence of a large excess of an 
RNA inhibitor that corresponds to the portion of an 
RNA substrate that forms base-pairs with the IGS. 
As in our previous selection procedure, ribozymes 



34 


Specialization of a DNA-cleaving Ribozyme 


that cleaved the DNA substrate became tagged by 
attachment of the 3' portion of the substrate to their 
3’ end, and thus were selectively converted to 
cDNA and amplified. On the other hand, ribozymes 
that preferentially bound the RNA inhibitor were 
prevented from binding and cleaving the DNA 
substrate, and thus could not be amplified. 

Employing this selection scheme, we subjected 
the pool of ribozymes obtained after 27 generations 
of in vitro evolution to an additional 36 generations, 
selecting for DNA cleavage activity in the presence 
of progressively increasing amounts of the RNA 
inhibitor. Based on kinetic analysis of the resulting 
individuals, the population responded as expected. 
The ratio of DNA to RNA-cleavage activity 
improved about 100- fold, with about a 50-fold 
increase in k ca i/K m for DNA cleavage and a slight 
decrease in k c * »/K m for RNA cleavage, measured 
for the population as a whole. A characteristic 
individual isolated from the evolved population 
exhibited kc at /K m values of 1.6 x 10 7 M' 5 min" 1 and 
6.9 x 10* M" 1 min" 1 for the DNA and RNA-cleavage 
reactions, respectively. This altered behavior, 
relative to the DNA-cleaving ribozymes obtained 
earlier, is largely due to two concerted mutations 
that arose within the catalytic core of the ribozyme 
over the course of the in vitro evolution process. 

Results 

Evolutionary Improvement in DNA 
substrate specificity 

Beginning with the population of DNA-cleaving 
ribozymes obtained after 27 generations of in vitro 
evolution (Tsang & Joyce, 1994), we carried out 36 
additional generations of the evolution process, 
applying a new selection strategy that simul- 
taneously selected for cleavage of a DNA substrate 
and against binding of an RNA inhibitor (Figure 2). 
The RNA inhibitor, corresponding to the 5' product 
of the ribozyme-catalyzed RNA-cleavage reaction, 
bound tightly to the IGS of the ribozyme 
(JK d < 1 nM; Tsang & Joyce, 1994), thus preventing 
binding of the DNA substrate. In this way, 
ribozymes that preferentially bound RNA over 
DNA were selectively disfavored. At the outset, the 
concentration of DNA substrate and RNA inhibitor 
was 0.2 pM and 0.1 pM, respectively. The concen- 
tration of RNA inhibitor was increased gradually 
over the course of in vitro evolution to 2 pM in order 
to maintain strong selection pressure against RNA 
binding. 

Selection against RNA binding also could have 
been accomplished using a full-length RNA 
substrate as the inhibitor. However, because the 
chemical step of RNA cleavage by the wild-type 
ribozyme is very fast and substrate binding is rate 
limiting (Herschlag & Cech, 1990), we felt that the 
same negative selection pressure could be exerted 
using either the RNA substrate or the 5' RNA 
product. Carrying out selection in the presence of 



generation 

Figure 3. Change in the catalytic efficiency of DNA and 
RNA-cleavage reactions for the evolving population of 
ribozymes. The k ca i/K m values for DNA cleavage (filled) 
and RNA cleavage (open) were determined under 
single-turnover conditions for the population as a whole, 
based on a least-squares fit to the equation: k obs - (k cal / 
K«n)[E], where [E] « K m and k& $ is the initial rate of the 
reaction (see Materials and Methods). PCR mutagenesis 
was performed twice at G36 and three times at G45 and 
G54 to increase the diversity of the evolving population, 
as described in Materials and Methods. 


the 5' RNA product provided additional negative 
pressure against ribozymes that catalyzed a ligation 
reaction resembling the reverse of the DNA-cleav- 
age reaction, in which the 5’ RNA product bound 
to a reacted DNA-cleaving ribozyme and attacked 
the ligation junction between the 3’ end of the 
ribozyme and the 3‘ portion of the DNA substrate, 
thereby releasing the tag that identifies the 
ribozyme for selection. 

In order to assess the impact of our selection 
strategy, we determined the catalytic efficiency of 
both DNA and RNA cleavage for the population of 
ribozymes obtained after the 27th, 36th, 45th, 54th, 
and 63rd generations of in vitro evolution (Figure 3). 
The preference for DNA over RNA, as determined 
by the relative k cM /K m values for the reaction with 
DNA versus RNA substrates, increased steadily, 
improving 100-fold between generations 27 and 63. 
There was an overall 50-fold increase in DNA- 
cleavage activity and a twofold decrease in 
RNA-cleavage activity. From generation 27 through 
36, there was a modest improvement in DNA- 
cleavage activity, but no significant change in 
RNA-cleavage activity. Based on sequence analysis 
(see below), we concluded that the population at 
that point in the evolution process lacked sufficient 
diversity for the desired phenotype to emerge. 
Accordingly, we subjected the population at 
generation 36 to two successive rounds of 
mutagenic PCR, introducing random mutations at 
an expected frequency of 1.4% per nucleotide 
position (Cad well & Joyce, 1992). Following this 
randomization procedure and an additional nine 
generations of in vitro evolution, RNA-cleavage 
activity declined significantly (Figure 3). Over the 
final 18 generations, from generations 45 through 
63, both DNA and RNA-cleavage activity rose 
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Table 1. Catalytic properties and sugar-specificity of DNA-cleaving ribozymes 






DNA 


RNA 


Ribozyme 

Mutations 

Positions 

271/312 

k (M e (min* 1 ) 

Kw (nM) 

kn it / Km 

(IQ 7 M _1 min' 1 ) 

k< u /K m * 
O0 7 M-’ min' 1 ) 

Specificity 

(DNA/RNAV 

Wild-type* 

0 

u/c 

2.4(4- 0.2) x 10~* 

6.0 ± 1.7 

4.0( + 1.6)xl0^ 

0.97 4- 0.19 

4.1 ( 4-2.6) x 10“* 

G27 48 b 

17 

U/A 

0.67 ± 0.02 

0.21 ± 0.03 

032 ± 0.05 

2.6 ± 0.05 

0.12 ± 0.02 

C45 8 

21* 

C/AA 

0.21 ± 0.01 

0.043 ± 0.007 

0.48 ± 0.1 

0.60 4- 0.02 

0.81 ± 0.2 

G6319 

Site-directed 

mutants: 

29' 

C/AA 

0.26 ± 0.01 

0.016 ± 0.002 

1.6 ± 0.2 

0.69 ± 0.02 

2.4 ± 0.3 

G63 19 

29 

C/A 8 

0.10 ± 0.00 

0.028 ± 0.003 

0.35 4- 0.05 

0.43 4- 0.02 

0.82 + 0.2 

G63 19 

28 

C/G 8 

0.021 ± 0.001 

0.0092 ± 0.002 

0.23 ± 0.05 

0.30 ± 0.07 

0.75 ± 0.4 

G6319 

28 

U/AA h 

1.5 ± 0.1 

0.33 4- 0.07 

0.47 ± 0.1 

1.2 ± 0.08 

0.40 ± 0.2 

G6319 

28 

U/A 8 

1.1 ± 0.1 

0.13 + 0.03 

0.94 ± 0.4 

1.4 ± 0.04 

0.67 ± 0.3 

G6319 

27 

U/G* 

0.44 4- 0.02 

0.042 4- 0.005 

1.1 ± 0.2 

1.3 ± 0.03 

0.82 ± 0.1 


* Data obtained previously (Beaudry & Joyce, 1992), modified slightly as a result of subsequent statistical analysis. 

b Data obtained previously (Tsang & Joyce, 1994). 

( Determined under single- turnover conditions ((El » [S]), as described in Materials and Methods. Data were fit to the equation: 
k» b* = k<M [E] /([El + K m ), which holds true when [S] « K m ; was obtained from the equation: y = l - e" k "* wl , where y = fraction 
substrate cleaved at time / for a range of [EJ spanning K m . 

d Determined under single-turnover conditions, as described in Materials and Methods. Data were fit to the equation: fcoi* « 
(Jr«t/K m ) [E], where k is the initial rate of the reaction at varying [E] when IE] « K m . 

* 44:G -♦ A, 51/52:CUAA insert, 86/87:UA - del, 94:A - U, 115/116:AG - UA, 165:C - U, 170:C - U, 188:G - A, 190/ 
191 :UG -* AU, 205 :U - C, 215:G - A, 239:U - A, 271:U - C, 289:U - A, 312:C - AA, 340:U - A, 364:C - U, and 366:G - A. 

f 44:G -> A, 51 :G - A, 51/52.CUAA insert, 86 /87:UA - del, 94:A - U, 115/116:AG - UA, 138:C - A, 165:C - U, 170:C - U, 
175:G -♦ A, 188:G -* A, 190/191:UG -» AU, 205:U - C, 215:G - A, 239:U - A, 271:U^C, 288/289:GU -♦ AA, 292:U -* A, 
312:G - AA, 334:A - G, 340:U - A, 347:A - G, 366:G - A, 389:U - C, and 390:A - C. 

* Mutations introduced by template mutagenesis (see Materials and Methods). 

h Mutations introduced by nested PCR. 

* Ratio of fccai/Km for DNA cleavage versus RNA cleavage. 


steadily. By generation 63, k^/K m values had 
reached 8 x l(r M' 1 min' 1 and 9 x 10 6 M' 1 min' 1 for 
the DNA and RNA-cleavage reactions, respectively, 
measured for the population as a whole. 

Comparison of individual ribozymes 

Individual ribozymes were isolated from the 
population at nine-generation intervals from gener- 
ation 27 through 63 using a shotgun cloning 
procedure (see Materials and Methods). Represen- 
tative individuals were chosen from generations 27, 
45, and 63 for formal kinetic analysis (Table 1). 
These ribozymes were representative of the 
population from which they were cloned because 
they contained the typical number and constella- 
tion of mutations that occurred in the population. 

The G63 19 ribozyme catalyzed the DNA- 
cleavage reaction with a k C9t /K m value of 
1.6 x 10 7 M' 1 min -1 , about fivefold more efficiently 
than the G27 48 ribozyme. This higher catalytic 
efficiency was due to a 13-fold lower K m that was 
partially offset by a 2.6-fold lower k cM . In contrast, 
the G63 ribozyme catalyzed RNA cleavage about 
fourfold less efficiently compared with the G27 
ribozyme. The preference for DNA versus RNA 
substrates, as determined by the ratio of k ca{ /K m 
values for the respective cleavage reactions, 
improved about 20-fold, from 0.12 for the G27 
ribozyme to 2.4 for the G63 ribozyme. The 
DNA-cleavage reaction proceeded to completion 
when carried out in the presence of an excess 
amount of either the G27 or G63 ribozyme. In 
contrast, the RNA-cleavage reaction did not 


proceed to completion, reaching about 40% cleaved 
for the G27 ribozyme and only about 15% cleaved 
for the G63 ribozyme (Figure 4). 

Evolutionary changes in genotype 

In order to determine which genotypic changes 
were responsible for the improved phenotype of 
the evolved population, sequence data were 
obtained from 25 randomly chosen individuals, 
isolated from the population every ninth generation 
from G27 through G63, as summarized diagram- 
matically in Figure 5. The mean number of 
mutations per cloned individual, relative to the 
wild-type, increased from 17 at G27 to 27 at G63, 
reflecting the continual accumulation of mutations 
that occurred during the in vitro evolution process 
(Figure 6(a)). However, the genetic diversity of the 
population, as measured by the Shannon entropy 
statistic (see Materials and Methods), declined 
steadily between GO and G36. In order to increase 
diversity, new mutations were introduced at an 
enhanced rate through successive mutagenic PCR 
procedures, performed at G36, G45 and G54. This 
resulted in increased genetic diversity between G45 
and G63 (Figure 6(b)). 

Prior to generation 18, mutations accumulated in 
both the phylogenetically conserved portions of the 
molecule that were randomized in the initial 
population and the peripheral regions that were not 
randomized (Figure 6(a)). After G18, however, 
mutations tended to accumulate in the periphery 
but not in the conserved core. Accordingly, the 
peripheral regions became more genetically diverse 
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(a) 




Figure 4. Comparison of the extent of the DNA and 
RNA-cleavage reactions catalyzed by the G27 48 and 
G63 19 ribozymes. Reactions were carried out in the 
presence of 23 (aM ribozyme and trace amounts (<1 nM) 
of either RNA or DNA substrate. (O) DNA cleavage by 
G27 48; (•) DNA cleavage by G63 19; (□) RNA cleavage 
by G27 48; (■) RNA cleavage by G63 19. (a) Time course 
over 120 minutes, (b) Expanded view of the reaction over 
the first two minutes. Data for RNA cleavage and DNA 
cleavage by G27 8 were fit to a biphasic exponential 
equation with a fast phase followed by a slow phase. 
Data for DNA cleavage by G63 19 were fit to a single 
exponential equation. All curve fits were based on the 
complete data set obtained over 120 minutes. 


while the core became more homogenous (Fig- 
ure 6(b)). By generation 63, the majority of 
mutations occurred in the periphery. 

Two major genotypic changes first appeared in 
the population at G36. The first, a CUAA insert 
located between nucleotide positions 51 and 52 
within the P2 stem (Figure 1), was present in >90% 
of the clones at G36, replacing the NGAA insert that 
had occurred with a frequency of 30% at G27. The 
second, a C -► U mutation at position 165 within 
P5c was present in >90% of the clones at G36, 
replacing a C -► A change at the adjacent position 
166 that had occurred in nearly all of the clones at 
G27. Between G36 and G45, five new mutations 
became dominant (>90% frequency) in the popu- 
lation. Three of these mutations, 271 :U -» C, 
289:U -> A, and 312:G -► AA, occurred within the 
catalytic core of the ribozyme and two, 340:U -► A 
and 364:C -+ U, occurred in the peripheral P9 
stem. Between G45 and G63, no new mutations 
rose to dominance in the population. Yet, the 


mean number of mutations per cloned individual 
continued to rise as a result of various low- 
frequency changes that occurred in the P9 and L8 
peripheral regions. 


Analysis of specific mutations 

Two of the mutations within the catalytic core, 
271:U -> C and 312:G -> AA, were chosen for more 
detailed analysis. Both of these changes were 
observed in all 25 clones that were obtained at G45 
and G54 and in 24 of 25 clones that were obtained 
at G63. No individual was found that contained 
only one of these two mutations, suggesting that 
they operate in concert to confer selective advan- 
tage to the ribozyme. The single individual at G63 
that lacked the two mutations contained instead the 
wild-type nucleotide U at position 271 and a G A 
change at position 312. Previous studies had shown 
that the 312:G -► A mutation occurred in a 
mutually exclusive fashion with regard to the 
concerted mutations 313;G -► Y and 31 4: A -+ G 
(Beaudry & Joyce, 1992; Tsang & Joyce, 1994). Over 
the course of evolution, 313:G -» Y and 314:A -* G 
nearly always occurred together, as their frequency 
rose, then fell, then rose again, and finally fell to 
extinction (Figure 7). The 312 mutation similarly 
rose and fell but in opposite phase relative to the 
313:G -> Y and 314:A -► G mutations. Prior to G36, 
the 312 mutation occurred as a G -♦ A change, but 
following G36 it occurred as a G -► AA change. 

Against the background of the G63 19 ribozyme, 
we reverted the 271 :U -* C mutation to the 
wild-type U and the 312:G -> A A mutation to either 
the wild-type G or a single A. All pairwise 
combinations of these mutations were prepared, 
and their effect on the RNA-catalyzed phosphoester 
transfer reaction involving either an RNA or DNA 
substrate was analyzed (Table 1). Reversion of both 
mutations to the wild-type state resulted in only a 
slight decrease in the efficiency of DNA cleavage, 
while reversion of either mutation alone led to 
a more substantial decrease (Figure 8). The 
271 :U -► C mutation, regardless of whether the 312 
position was G, A or AA, resulted in lower 
RNA-cleavage activity. The 271 :U -► C mutation 
also resulted in a lower k Cl , t value for DNA cleavage, 
but this was compensated by a much lower K m 
value. The 312:G -► A A change resulted in a sig- 
nificantly higher k cal value for DNA cleavage, 
regardless of whether the 271 position was U or C. 
However, for 271 :U, but not for 271 :C, this was 
accompanied by a substantial increase in K m value. 
In summary, the 271 :U C and 312:G -*• A A 
mutations were synergistic: together they provided 
the greatest improvement in DNA-cleavage activity 
while diminishing RNA-cleavage activity. 

This mutational analysis clearly demonstrates the 
selective advantage of the concerted nucleotide 
mutations at positions 271 and 312. These two 
changes, however, do not explain all of the 
improvement in DNA versus RNA-cleavage activity 
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Figure 5. Sites at which mutations occurred over the course of evolution, superimposed on the secondary structure of the Tetrahymcna ribozyme. Box height corresponds 
to the frequency of mutations (%) at each nucleotide position, on the basis of 50 subclones sequenced at G27, and 25 subclones sequenced at G36, G45, G54, and G63. 
Non-mutable primer binding sites are shaded; substrate is shown in black. Commonly occurring mutations ( >50% frequency) are labeled, (a) G27 (Tsang & Joyce, 1994); 
(b) G27, normalized such that mutations occurring with > 50% frequency were assigned a value of 0% in this and subsequent panels; (c) G36; (d) G45; (e) G54; (f) G63. 
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that occurred between G27 and G63. Even when 
both mutations were reverted to the G27 state 
(271 :U; 312: A) on the G63 19 background, the ratio 
of DNA to RNA-cleavage activity was still 
significantly greater than that of the G27 48 
ribozyme. Therefore, other mutations that arose 
between G27 and G63 must contribute to the 
enhanced substrate specificity of the evolved 
ribozymes, such as those that occurred in the P9 
and L8 regions. 

Discussion 

Preference for DNA versus RNA substrates 

We have described a novel in vitro evolution 
strategy that simultaneously selects for cleavage of 
DNA and against binding of RNA, leading to 
specialization of the DNA-cleaving activity of a 
group I ribozyme. Beginning with a heterogeneous 
population of ribozymes that cleave both RNA and 
DNA substrates, we applied this strategy for 36 
generations of in vitro evolution and obtained a 
50-fold increase in DNA-cleavage activity and a 
twofold decrease in RNA-cleavage activity, corre- 
sponding to an overall increase of 100-fold in the 
preference for DNA over RNA substrates. From the 
evolved population at G63, we isolated a ribozyme 
that cleaves DNA more efficiently than it cleaves 
RNA, with a catalytic efficiency, k ca t/K m , of 
2 x 10 7 M' 1 min' 1 and 7 x 10 6 M' 1 min' 1 for the DNA 
and RNA-cleavage reactions, respectively. Site-di- 
rected mutagenesis experiments confirm the selec- 
tive advantage of two covarying mutations at 
nucleotide positions 271 and 312 within the 
catalytic core of the ribozyme that are largely 
responsible for this modified behavior. 

Our selection regime required that the ribozymes 
cleave a DNA substrate in the presence of an RNA 
inhibitor that corresponds to the 5’ product 
generated by the RNA-catalyzed RNA-cleavage 
reaction. Ribozymes with enhanced affinity for the 
DNA substrate were more likely to avoid inhibition 
by the RNA product. Under pseudo-first-order, 
single-turnover conditions (see Materials and 
Methods), substrate binding is reflected by K m , the 
concentration of ribozyme that is required to obtain 
the half-maximal first-order rate of reaction. The 
G63 19 ribozyme catalyzed the DNA-cleavage 
reaction with a K m value that is 13-fold lower than 
that of the G27 ribozyme. Thus, the G63 ribozyme's 
enhanced preference for DNA versus RNA sub- 
strates compared to the G27 ribozyme may be due, 
at least in part, to the G63 ribozyme's greater 
affinity for the DNA substrate. 

Ribozymes that bound the RNA product more 
slowly also would have had an increased likelihood 
of escaping inhibition by the RNA product under 
our selection conditions, assuming that the rate of 
DNA substrate binding did not decrease as well. 
Changes in RNA binding rates may explain the 
differing maximal extent of RNA cleavage by the 


(a) 




Figure 6. Genotypic change over 63 generations of 
in vitro evolution, (a) Mean number of mutations per 
individual, pertaining to the entire ribozyme (#) or the 
140 nucleotides within the catalytic core that were 
partially randomized in the initial population (Ok Data 
for G9 are from Beaudry & Joyce (1992); data for G18 and 
G27 are from Tsang & Joyce (1994). (b) Shannon diversity 
of randomized core (O) and peripheral regions (□), 
calculated as described in Materials and Methods. 


G27 and G63 ribozymes: G27 48 cleaves RNA to a 
maximal extent of ^-40%, while G6319 reaches 
only -15%. For the wild-type Tetrdhymena ri- 
bozyme, incomplete cleavage of the RNA substrate 
occurs because an internal equilibrium is estab- 
lished between the forward and reverse reactions 
due to slow dissociation of the product (Mei & 
Herschlag, 1996). A similar situation may exist for 
the G27 and G63 ribozymes. If the rate of the 
reverse reaction is similar for these two ribozymes, 
then the lesser extent of RNA cleavage by G63 19 
compared with G27 48 is likely due to the slower 
rate of forward reaction for the former. 

Genetic basis for the evolved phenotype 

Site-directed mutagenesis experiments indicate 
that covarying mutations at positions 271 and 312 
in the G63 19 ribozyme contribute to its preference 
for DNA versus RNA substrates (Table 1). Both of 
these mutations disrupt phylogenetically conserved 
base-pairing interactions within the catalytic core 
(see Figure 1). The 271 :U — ► C mutation disrupts 
what is thought to be a reverse Hoogsteen pair 
between 103 A and 271 U, located at the end of the 
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313/314 312 

GA-*CQ G~*AA 



Figure 7. Frequency of mutually exclusive mutations at 
positions 312 (•) and 313/314 (O) over the course of 
in vitro evolution. Data obtained at discrete times were fit 
to a spline curve. The 312:G -* A mutation was not 
detected between G27 and G36, but reappeared 
beginning at G45 as 312:G -+ AA. The 313/314:GA -+ UG 
mutation was replaced with 313/314:GA -► CG begin- 
ning at G27, but was not detected at G45 and thereafter. 


P3 stem (Kim & Cech, 1987; Green el <?/., 1990). The 
312:G AA mutation disrupts the 262C-312G 
Watson-Crick pair, which lies within the P7 stem 
adjacent to the 264G-311C base-pair that forms a 
base-triple with the bound guanosine nucleophile 
(Michel et al., 1989). By disrupting the local 
structure near the G-binding site, the 271 and 312 
mutations may alter the orientation of the attacking 
3'-terminal guanosine relative to the target phos- 
phate, differentially affecting the DNA versus 
RNA-cleavage reactions. Although covarying mu- 
tations are often observed between two interacting 
components of an enzyme, the nucleotides at 
positions 271 and 312 are unlikely to interact 
directly because they are located at opposite ends 
of the P7 stem. Their codependent behavior more 
likely results from indirect interactions. Other 
mutations that contribute to the altered phenotype, 
including those in the L8 and P9 regions, occur 
outside the catalytic core of the ribozyme, 
suggesting that they too exert their influence via 
indirect or long range effects. 

Residues in the J8/7 region of the wild-type 
Tetrahymena ribozyme have been shown to interact 
directly with 2'-hydroxyl groups of the substrate 
strand of the PI duplex (Pyle et al., 1992). Thus, we 
expected that selection against RNA binding would 
lead to mutations in this region. Yet, J8/7 proved 
to be highly resistant to mutation over the course 
of in vitro evolution. Not a single individual among 
the 100 that were isolated from generations 36, 45, 
54, and 63 had a mutation in this region. 
Furthermore, when we intentionally mutated 
nucleotide position 302, which in the wild-type 
ribozyme interacts with a particular 2'-hydroxyl of 
the RNA substrate, we observed a near-complete 
loss of DNA-cleavage activity (data not shown). 
Therefore, despite the absence of 2'-hydroxyl 



position 271/312 

Figure 8. Co-dependence of the 271 :U — ► C and 
312:G -*• AA mutations with regard to enhanced prefer- 
ence for DNA versus RNA substrates. The k <a t/K m values 
for the DNA (filled) and RNA (open) cleavage reactions 
were determined for ribozymes containing neither, 
either, or both of the 271 and 312 mutations in the context 
of the G63 19 ribozyme. 


groups in the DNA substrate, specific nucleotides 
within J8/7 appear to be critical for DNA-cleavage 
activity. 

Although the G63 19 ribozyme contains 29 
mutations relative to the wild-type sequence 
and cleaves RNA fourfold less efficiently than the 
G27 48 ribozyme, it still cleaves RNA quite readily, 
with a kcnt/Km value of 6.9 x 10 6 M" 1 min' 1 , which is 
similar to the wild-type level of 9.7 x 10 6 M' 1 min -1 . 
The ribozyme's ability to tolerate a large number of 
mutations suggests that it is overspecified with 
respect to RNA-cleavage activity. More dramatic 
sequence or structural alterations may be necessary 
in order to achieve a substantial reduction in this 
activity. However, elimination of RNA-cleavage 
activity without significant reduction of DNA- 
cleavage activity is likely to be difficult to achieve 
for several reasons. Firstly, the ribose at the 
cleavage site is inherently more reactive than 
deoxy ribose, contributing, in the case of the 
wild-type Tetrahymena ribozyme, ^10 3 -fold 
towards stabilization of the transition state com- 
pared with deoxyribose (Herschlag et al., 1993). 
Secondly, RNA/RNA homoduplexes are more 
stable than DNA/ RNA heteroduplexes, resulting 
in inherently greater binding affinity for RNA 
compared to DNA substrates (Roberts & Crothers, 
1993). Thirdly, both RNA/RNA and RNA/DNA 
duplexes tend to adopt A-form helices (Wang et aL, 
1980; Sanger, 1984; Egli et al., 1992), making it 
difficult for the ribozyme to distinguish them, 
especially if the ribozyme is to maintain generality 
with respect to the sequence of the ribozyme /sub- 
strate pairing. Fourthly, steric conflicts involving 
2’-hydroxyl groups, which might reduce RNA- 
cleavage activity, may not be substantial due to 
the small size of these groups. Considering the 
inherent advantages of a ribose compared to 
deoxyribose-containing substrate, it is remarkable 
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that the latter is cleaved more efficiently by the 
evolved ribozyme. 

Evolutionary search for fitness optima 

In intro evolution experiments provide an 
opportunity to examine the genotypic and pheno- 
typic responses of an evolving population to 
controlled selection pressures. Our test-tube popu- 
lation of ribozymes displayed some of the 
tendencies of other artificially and naturally-evolv- 
ing populations, including the occurrence of 
synergistic and mutually exclusive traits. For 
example, over the course of 63 generations of 
evolution, the DNA-cleaving ribozymes generally 
possessed either a mutation at position 312 or a 
pair of concerted mutations at positions 313 and 
314. These two alternatives were mutually exclu- 
sive, alternately rising and falling in abundance 
(Figure 7). Lehman & Joyce (1993) observed a 
similar response in an evolving population of 
Ca 2+ -dependent RNA-cleaving ribozymes. Such 
complex behavior demonstrates that artificially- 
evolving populations, like natural populations, 
exhibit an adaptive response that fluctuates as 
variants of variants are generated over the course 
of evolution. 

Despite increasing diversity in the evolving 
population between generation 45 and 63, pheno- 
typic improvements were modest and genotypic 
changes occurred predominantly in the peripheral 
regions of the ribozyme. Few mutations were 
acquired in the critical core regions, where diversity 
remained low (Figure 6(a) and (b)). These results 
suggest that the evolved population may have 
become trapped in a local fitness optimum. 
Achieving a higher level of fitness might require 
traversing less advantageous intermediates, for 
example, by exploring a much wider distribution of 
variants or by restarting the evolution process at an 
earlier generation. Alternatively, the continued 
accumulation of neutral mutations, in peripheral 
regions, for example, might enable the population 
to bypass less favorable intermediates and reach a 
state of improved fitness through a global change 
in the secondary structure of the ribozyme 
(Schuster et al ., 1994; Schuster, 1995). 

In this study, we demonstrated that by selecting 
for DNA cleavage and against RNA binding, one 
can enhance the specificity of a group I ribozyme 
for DNA versus RNA substrates. Greater specificity 
might be obtained by carrying out additional 
rounds of in vitro evolution under increasingly 
stringent selection conditions, perhaps employing 
higher concentrations of RNA inhibitor or shorter 
incubation times. Alternatively, increased prefer- 
ence for DNA versus RNA substrates might be 
achieved by focusing the selection pressure 
exclusively on the DNA-cleavage reaction, seeking 
maximal activity with a DNA substrate and 
expecting, as some have argued (Eaton et al ., 1995), 
that improved specificity will follow automatically. 
Because the evolved DNA-cleaving ribozymes are 


already efficient catalysts, adopting this strategy 
will require the introduction of extremely stringent 
selection constraints that may be technically more 
difficult to implement. Of course, whether further 
improvement is even possible depends ultimately 
on the functional flexibility of the group I catalytic 
motif and the catalytic potential of RNA. 

Materials and Methods 

Materials 

Unlabeled nucleoside triphosphates (NTPs) and de- 
oxynucleoside triphosphates (dNTPs) were purchased 
from Pharmacia, and dideoxynucleoside triphosphates 
(ddNTPs) were from U.S. Biochemical. [a- 32 P]GTP, 
[y- 32 P]ATP, and [ 3 H]UTP were from ICN Radiochemicals. 
Synthetic oligodeoxynucleotides were obtained from 
Operon Technologies and purified by polyacrylamide gel 
electrophoresis and subsequent chromatography on 
either Sephadex G-25 or DuPont Nensorb 20. Synthetic 
oligoribonucleotides were obtained from Oligos Etc. and 
purified by gel electrophoresis and subsequent chroma- 
tography on DuPont Nensorb. Bacteriophage T4 polynu- 
cleotide kinase was from New England Biolabs, calf 
intestine phosphatase from Boehringer, avian myeloblas- 
tosis virus (AMV) reverse transcriptase from Life 
Sciences, and Moloney murine leukemia virus (MoMLV) 
reverse transcriptase, Sequenase 2.0 (modified bacterio- 
phage T7 DNA polymerase), T4 DNA ligase, and T4 
DNA polymerase from U.S. Biochemical. T7 RNA 
polymerase was prepared as described (Davanloo et al., 
1984) and purified according to a procedure originally 
developed for SP6 RNA polymerase (Butler & Chamber- 
lin, 1982). Streptavidin immobilized on agarose beads 
was obtained from Sigma and streptavidin-containing 
pipette tips (Affinitips) were from Genosys. 

In vitro evolution protocol 

Starting with the G27 population, in vitro evolution 
was carried out as described (Tsang & Joyce, 1994), with 
the following modifications. (1) The DNA-cleavage 
reaction was performed in the presence of an RNA 
inhibitor, consisting of the 5’-product portion of the RNA 
substrate, 5-GGCCCUCU-3'. Generally, 0.1 to 2% of the 
population underwent the target reaction and was 
selectively amplified. The concentration of inhibitor was 
increased when the reacted portion of the population 
surpassed 2%, in order to maintain a strong selective 
advantage for the more substrate-specific DNA-cleaving 
ribozymes. Over the 36 generations the inhibitor 
concentration was increased from 0.1 to 2 |iM. (2) At G37, 
the ribozyme concentration was increased from 0.1 nM to 
1 nM, and the incubation time was increased from five 
minutes to one hour in order to give ribozymes with 
lower activity but higher specificity the opportunity to 
react with the DNA substrate. (3) In order to prevent 
isothermal amplification artifacts, reacted ribozymes 
were purified away from non-reacted molecules using 
either of two methods. From generation 28 through 36, 
the ribozyme population was challenged with a 
3'biotinylated DNA substrate so that reacted molecules 
became tagged with the biotinylated portion of the 
substrate. Active molecules were isolated by binding to 
streptavidin immobilized on agarose beads. Non-specific 
binders were removed by extensive washing with a 
solution containing 1 M NaCl, 50 mM Tris (pH 7.5), 5 mM 




Specialization of a DNA-cleaving Ribozyme 


41 


EDTA. The agarose-bound molecules were then am- 
plified in an isothermal reaction, as described (Tsang & 
Joyce, 1994). The agarose beads did not significantly 
affect the efficiency of amplification. From G37 through 
G63, reacted ribozymes were separated from non-reacted 
ones by denaturing polyacrylamide gel electrophoresis, 
then eluted from the gel and purified by chromatography 
on DuPont Nensorb. (4) New mutations were introduced 
into the population at generations 36, 45, and 54 through 
mutagenic PCR (Cadwell & Joyce, 1992; Cadwell & Joyce, 
1994). The mutagenic PCR was seeded with cDNA 
obtained by reverse transcription of RNA that had been 
purified by gel electrophoresis. The same primers were 
used in both the standard and mutagenic PCR: 
5-CG AGTACTCC A A AACT A ATC-3' (primer lb) and 
5-CTGCAGAATTCTAATACGACTCACTATAGGAGG- 
GAAAAGTTATCAGGC-3* (primer 2). At G36, two con- 
secutive rounds of mutagenic PCR were carried out, pro- 
viding an expected mutation rate of 1.4% per nucleotide 
position; at G45 and G54, three consecutive rounds were 
carried out, providing an expected mutation rate of 2% 
per position. In each case, the cDNA strand was isolated 
and purified by polyacrylamide gel electrophoresis prior 
to seeding the next round of mutagenic PCR. (5) Our 
standard PCR protocol was altered slightly between G45 
and G63, reducing the number of temperature cycles 
from 30 to 15, and carrying out two consecutive 15-cycle 
amplification reactions. This was done in order to prevent 
PCR artifacts and to reduce the amount of carryover 
material from one generation to the next. 

Kinetic analysis 

All kinetic analyses was carried out under single- 
turnover conditions with ribozyme (E) in large excess 
over substrate (S), and the concentration of S at least 
tenfold below K m . Ribozymes were preincubated for 15 
minutes at 37°C in a 17 pi solution containing 10 mM 
MgCU and 30 mM 4-(2-hydroxyethyl)piperazine-l-pro- 
panesulfonic acid (EPPS) (pH 7.5). A total of 3 pi of 
substrate in the same solution was added to initiate the 
reaction, and 3 pi aliquots were removed at the appro- 
priate times and quenched by the addition of one to two 
volumes of an ice-cold mixture containing 8 M urea, 
50 mM to 100 mM EDTA, 0.05% (w/v) xylene cyanol, 
0.05% (w/v) bromophenol blue, 10% (w/v) SDS, 9 mM 
Tris-borate (pH 8.3), and 20% (w/v) sucrose. The k cM /K m 
values were calculated from the slope of the best fit line 
obtained by plotting k oin versus [E] for at least four 
different subsaturating concentrations of E. The k cM and 
K m values were determined individually for the DNA- 
cleavage reaction by plotting k 0 over a range of ribo- 
zyme concentrations that spanned K m . Data were fit to 
the equation: k<*» = k c »t /([E] 4- K m ), which holds true when 
[S] « K m , assuming that binding of the ribozyme and 
substrate has reached pre-equilibrium. The standard 
errors reported for k cM and K m values were determined 
based on the fit to this equation by non-linear regression 
using the Marquardt-Levenberg algorithm (Sigma Plot). 
The k CM value for RNA cleavage was considerably faster 
than 1 min" 1 and thus not measurable using manual 
pipetting methods. Determination of kinetic parameters 
on different days and with different preparations of 
material generally varied by less than 20%. DNA and 
RNA-cleavage assays were carried out using the same 
preparation of ribozyme transcript and generally in a 
side-by-side fashion to obtain the most accurate 
measurement of the ratio of DNA to RNA-cleavage 
activity. 


Preparation of ribozymes and substrates 

Individual ribozymes were isolated by shotgun cloning 
using the Invitrogen TA Cloning Kit. Sequencing of 
clones and preparation of individual ribozymes were 
carried out as described (Tsang & Joyce, 1994). 
Preparation of DNA and RNA substrates and procedures 
for cDNA synthesis, PCR, and in vitro transcription 
reactions all were performed as described. 

Site-directed mutagenesis 

Mutations were made at specific positions in the 
ribozyme using a template-directed mutagenesis pro- 
cedure (Joyce & Inoue, 1982). The non-template strand 
was prepared by amplifying the ribozvme-encoding gene 
by PCR using a 5'-biotinylated primer for the template 
strand and a standard primer for the non-template 
strand. The PCR products were purified by phenol 
extraction and ethanol precipitation, then bound to a 
streptavidin-containing matrix (Affinitip) and washed 
with 0.2 N NaOH to denature and elute the non-biotiny- 
lated, non-template strand. The eluate was neutralized 
with HC1, and the non-template strand was concentrated 
and purified by ethanol purification. The 5'-phosphory- 
lated oligos that introduce the desired mutations 
(5 , -CCCCGACCGACATTTAGTCTGTGAACTGC-3‘ and 
5' - CTCCCATTAAGGAGAGGTCVG ACTATATCTTAT - 
GAGAAG-3’, Y = C or T) together with an oligo that 
defines the 3’ end of the template (primer lb) were 
annealed to the non-template strand, and the gaps were 
filled with T4 DNA polymerase and ligated with T4 DNA 
ligase. RNA was transcribed from the reconstructed 
double-stranded DNA, reverse transcribed to cDNA, 
amplified by PCR, subcloned, and sequenced. Individual 
ribozymes that contained the appropriate mutation were 
prepared in the usual manner. In the case of the 
271:C-+U mutation, site-directed mutagenesis was 
carried out using a nested PCR strategy. First, PCR was 
performed using the mutagenic oligo 5'^CCCCGACCG A- 
CATTTAGTCTGTGAACTGC-3 , and primer 2. The 
— 250nt product was purified using a Microspin S-30O 
HR column (Pharmacia) and ethanol precipitated. This 
product was used as primer in a subsequent PCR reaction 
with primer lb. The full-length (^400nt) PCR product 
was purified in a low-melting point agarose gel and 
recovered by phenol extraction and ethanol precipitation, 
and the mutant ribozyme was prepared in the usual 
manner. 


Diversity estimation 

Diversity in a population of genotypes was estimated 
by averaging, over all 353 variable nucleotide positions, 
the diversity at each position using the formula: 
H = Ip, (In pi ), where pi is the frequency of a particular 
nucleotide base (A, C, G, or U). Diversity is maximal 
when A, C, G, and U occur with equal frequency. 
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